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Abstract 

A particular choice of the time function in the recently presented spherical solution by 
Dadhich [1] leads to a singularity free cosmological model which oscillates between 
two regular states. The energy-stress tensor involves anisotropic pressure and a 
heat flux term but is consistent with the usual energy conditions (strong, weak and 
dominant). By choosing the parameters suitably one can make the model consistent 
with observational data. An interesting feature of the model is that it involves blue 
shifts as in the quasi steady state model [2] but without violating general relativity. 
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Following the discovery of non-singular cylindrically symmetric perfect fluid exact cosmo- 
logical solution of the Einstein equation by Senovilla [3], some spherically symmetric non- 
singular models have been presented by Dadhich et.al [1,4]. These models have an energy- 
stress tensor with anisotropic pressure and heat flux but obeying the strong, weak and 
dominant energy conditions. The metric has a time function which can be arbitrarily cho- 
sen subject to the constraint of non-singularity and the energy conditions. It turns out that 
there exist different such choices which will be discussed in a detailed paper separately. In 
this letter we shall confine to the choice that gives an oscillatory behaviour of the universe 
without any singularity. 

The authors are not aware of any oscillatory singularity free model in classical general 
relativity (GR) while there are some oscillatory behaviour models proposed in the recent 
formulation of of quasi steady state cosmology (QSSC) [2]. As the QSSC models predict the 
possibility of blue shift, our model would also admit that possibility. Thus should obser- 
vations in future reveal blue shifts, it may simply indicate that the matter in the uinverse 
is not perfect fluid and one need not bring in the ideas of non-conservation as in QSSC 
contradicting GR. 

Our oscillatory model is described by the metric [1], 

ds * = (r 2 + p)dt 2 _ 2r2+ ^ r 2 _ r 2 (d 02 + sin 20 rf 2) (1) 

r z + P 

where P = P(t) which can be chosen arbitrarily. The choice P(t) = a 2 + b 2 cosuot with 
a 2 > b 2 will render oscillatory behaviour to the model without encountering divergence of 
any kinematical and physical parameters. In Ref. [1] the choice made was P(t) = a 2 + b 2 t 2 , 
which of course did not give oscillatory behaviour. 

The energy-stress tensor for imperfect fluid is given by [5], 

T ik = {p + p)uiU k - pg ik + Ap[dc k + ^{g ik - UiU k )} + 2qc {i c k) (2) 
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where Ui and q are respectively unit timelike and spacelike vectors, p energy density and 
p isotropic fluid pressure, Ap pressure anisotropy and the term involving q represents heat 
flux. 

We employ the comoving coordinates to write U{ = y/gooS® and take q = y^guSj. The 
kinematic parameters; expansion, shear and acceleration for the metric (1) read as follows: 

—Pr 2 2 r 

6 = 2(2r 2 + P)(r 2 + P)3/2' ^ = 3^' Ur = ~r^TP' (3) 

Now applying the Einstein equation, we obtain 
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(2r 2 + P) 3 /2( r 2 + p)- 

The pressure anisotropy Ap = p r — p± is given by 

(9r 2 + 5P)P 2 
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Now we choose for the time function 



(7) 
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P(t) = a 2 + b 2 cosut, a 2 > b 2 . 



(9) 



This lends oscillatory behaviour to the model which oscillates between the two regular states. 
The oscillation period is t = 2-k/uj, density is maximum at t — (2n + \)n /uj and it is min- 
imum at t — 2nn/uj for an integer n. The model could have as low and as high density 
as one pleases by choosing large values for the parameters a and b with the former being 
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as close to the latter from the above. The solution involves three parameters a, b and oo. 
Besides, when comparing with observational data, one needs a specification of the locale of 
observation, i.e. the time to of observation and r the radial coordinate of the observer. Since 
there is abundance of free parameters, it is therefore possible to coast the model as close to 
the observations as one pleases by suitale choice. 

Above all the most interesting feature of the model is oscillatory behaviour indicating like 
the steady state cosmology no beginning and no end. As in QSSC [2], this model would 
also predict blue shifts, should they be observed in future, one need not necessarily have to 
invoke non-conservation of energy but instead an imperfect fluid distribution could as well 
do without violating GR and the usual energy conditions. This is quite remarkable and 
interesting feature of our model. 

The pressure anisotropy and heat flux fall off as r~ 4 and they vanish at the centre r = 0. It 
is obvious that expansion and heat flux have similar behaviour (note that in Ref. [1], there 
was a sign error which indicated opposite behaviour), vanishing at out = 0,n and attaining 
maximum at ut = 7r/2,37r/2. That is in expanding phase heat flows out while the reverse 
happens for contracting phase. The pressure anisotropy could like density be made as small 
as one pleases and it changes sign at t — (2n+ l)/2nu. Acceleration vanishes at both ends, 
r — > 0, oo and is finite for all t at a given r. 

For the oscillatory choice (9), it can be easily checked that p > p r ,P± > 0, (p+p r ) 2 — 4q 2 > 
and p — p r — 2p± + \(p + p r ) 2 — Aq 2 } 1 / 2 > always. This ensures that all the (weak, strong 
and dominant) energy conditions are satisfied. All the physical and kinematic parameters 
always remain regular and finite.. The metric is simple enough to see that it is causally 
stable and geodesically complete. 

The present model gives a picture somewhat different from non-oscillating singularity-free 
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models [1,3,4,6]. There the universe has a state of infinite dilution both in the infinite past 
and future and in between there is a state of maximum contraction. In our case P(t) never 
becomes arbitrarily large and so there is no infinite dilution but periodically the physical 
variables oscillate between finite maxima and minima. However as r — > oo, all the physical 
variables tend to vanish as in other non-singular models. This behaviour is demanded by 
the theorem [7] that in non-singular models, the space average of all the kinematic scalars 
and physical parameters must vanish (for a weaker theorem on the vanishing of space time 
averages [8]. 

We have thus shown that it is possible to have a truely oscillatory singularity free spherical 
model within GR without violating conservation of energy and its usual conditions. Such a 
model would be consistent with blue shifts, should they be discovered in future. The model 
has a number of free parameters which can be suitably chosen to coast it arbitrarily close to 
the observations. One may not quite relish the introduction of an imperfect fluid but then we 
are not violating any stringent physical requirement. Abondoning the cosmological principle 
of isotropy and homogeneity requires hardly any apology but the introduction of a specially 
favoured centre of the universe may require some justification. However one cannot get the 
solution of the Einstein field equation without introducing some symmetry assumption. In 
any case the solution broadens our horizon about the potentialities of relativistic models. 
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